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^ . ABSTRACT 

' Far-UV spectroscopy from the FU SE satellite is analysed to uniquely probe spatial 

structure and clumping in the fast wind of the central star of the H-rich planetary 
nebula NGC 6543 (HD 164963). Time-series data of the unsaturated Pv AA1118, 
(~| \ 1128 resonance line P Cygni profiles provide a very sensitive diagnostic of variable 

^ |. wind conditions in the outflow. We report on the discovery of episodic and recurrent 

optical depth enhancements in the Pv absorption troughs, with some evidence for 
a 0.17-day modulation time-scale. SEI line-synthesis modelling is used to derive 
^ physical properties, including the optical depth evolution of individual 'events'. The 

' characteristics of these features are essentially identical to the 'discrete absorption 

components' (DACs) commonly seen in the UV lines of massive OB stars. We have 
also employed the unified model atmosphere code CMFGEN to explore spectroscopic 
signatures of clumping, and report in particular on the clear sensitivity of the Pv lines 
to the clump volume filling factor. The results presented here have implications for 
the downward revision of mass-loss rates in PN central stars. We conclude that the 
temporal structures seen in the Pv lines of NGC 6543 likely have a physical origin 
. that is similar to that operating in massive, luminous stars, and may be related to 

' near-surface perturbations caused by stellar pulsation and/or magnetic fields. 

' Key words: stars: outflows — stars: evolution — stars: individual: NGC 6543. 
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1 INTRODUCTION in this paper is to exploit new Far-Ultraviolet Spectroscopic 

Explorer {FUSE) satellite time-series spectroscopy of the 
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^ As the immediate precursors of white dwarfs, the central 1 i r ^t^^ r-r An ±. ■ i ^ ■ j_- 

, V , central star oi JNGC 6543 to uniquely examme the time- 
■ - ' stars or planetary nebulae (CSPNs) provide fundamental . , , , . . ,. r .. r , • i 

, , . variable characteristics ot its last wind, 

tests m key areas of stellar astrophysics, including evolution- 

, , T rnri J 1 i I, J 1 1 J An improved understandine of variability and clumping 

ary theory, non-Lib model atmospheres and plasma hydro- ^ ° j f a 

J . ». -i jijji-i in the fast winds of CSPN is important for several reasons: 

dynamics. An important challenge is to understand details ^ 

r n u iu 1 i u 1 i Temporally variable spatial structure implies a clumped and 

OI the mass-loss process through the planetary nebula stage. j f f f 

o -c 11 ii, £ i ■ J f noTjAT -J u f iu at least locally non-spherical fast wind, which can modify 

bpecmcally, the tast winds ot CbPN provide a probe ot the .j f > j 

, , J , 1 ui ii- f iu how the supersonic outflow interacts with the nebular mate- 

current mass-loss, and represent a valuable setting tor the ^ 

, A r J- J- J u-i-i i- ui i rial (e.g. Falle et al. 2002; Steffen & Lopez 2003). Substantial 

study ot radiative and mechanical interactions between stars ^ ° ' f / 

J . , . . , wind clumping can also significantly complicate determina- 

and their environments. fa o j t- 

o, ,1 r, 1 I, Til T i M- 1 TTii ' 1 i tlous of thc mass-loss rates from the central stars. For ex- 

Shortly after the launch ot the international Ultraviolet 

Explorer (lUE) satellite, it became clear that fast winds ^"^Pl^' '^''^''^ ^^^^l^^^^ °^ °Pti^^l ^'^d UV Unes in luminous, 

r i r /^onivT / tt i n^n n ■ massive O-type stars have provided evidence for extreme 

are a common feature of CSPN (e.g. iieap, 1979; Permotto, ^ 

rA i- o 1 T1 i- -inonN rni n/-i -TTir wiud clumpiug factors, with implications that the mass-loss 

Cerruti-Sola, Benvenuti, 1982). The P Cygni UV resonance ^ ' ^ 

ciu u jjji J • -J rates of OB stars may have previously been over-estimated 

line profiles have been decoded to record maximum wind j ^ j 

, ... r , r innn 1 -1 J i- j- Til by factors of between 3 to 10 (e.g. Hillier, Lanz & Heap 

velocities ot up to a few 1000 km s and estimates of the ■' \ o > 

, , , , in-ii -t in-6 i\/r -1 / 2003; Massa et al. 2003; Markova et al. 2004; Bouret, Lanz 

mass-loss rate range between ~ 10 to 10 JVl© yr (see ' ' ' ' 

e.g. Perinotto, 1989 and references within). Our objective ^ Sillier 2005; FuUerton, Massa & Prinja 2006). 

Recent Chandra and XMM— Newton observations 

have established the presence of diffuse X-ray emission from 
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(AWF); awxb@star.ucl.ac.uk (AWB) of the X-rays from the central star vicinity remains uncer- 



2 R.K. Prinja et al. 



Table 1. NGC 6543 central star parameters. 



Parameter 


Value 




Reference 


Sp. type 


Of-WR (H- 


rich) 


Mendez et al. (1990) 


Luminosity 


5200 L0 




de Koter et al. (1996) 


Tefl 


63000 K 




Georgiev et al. (2006) 


Radius 


0.6 R0 




Georgiev et al. (2006) 


Distance 


1001 ± 269 


pc 


Reed et al. (1999) 


Mass-loss rate 


~ 1x10"^ 


M0/yr 


de Koter et al. (1996); 








Georgiev et al. (2006) 


Terminal velocity 


1400 km s" 


-1 


This study. 


Wind flushing time 


45 mins 




This study. 




tain, and one possible origin would be the presence of shock- 
heated gas arising from instabilities in a variable fast wind 
(e.g. Akashi, Soker & Behar, 2006). 

The present study is also motivated by a need to un- 
derstand the relation between optical spectroscopic (radial 
velocity) and photometric changes evident in sca'ctmI CSPN. 
The frequently imaged non-spherical and bipolar morpholo- 
gies of the planetary nebulae (e.g. Zuckerman & AUer, 1986) 
raise the possibility that the planetary nebulae may have 
been ejected via common envelope interaction in close bi- 
nary central stars (e.g. Sandquist et al. 1998). Surveys of 
the PN binary population are currently frustrated by the 
physical interpretation of the observed radial velocity and 
photometric variations, which might at least in part be due 
to changes in the fast wind (e.g. Handler, 2003; De Marco 
et al. 2004). 

There is a clear need therefore to investigate further the 
detailed temporal variability of the fast winds in CSPN. To 
date, studies of variability in these outflows have primarily 
relied on (limited) multiple UV spectra obtained with the 
lUE satellite (e.g. Patriarchi & Perinotto, 1997). Gener- 
ally however the detection of wind line variability with lUE 
(or HST) was extremely difficult since the only UV reso- 
nance lines accessible are strongly saturated, thus masking 
all changes except at the extreme violet edges of the line pro- 
files. In addition, the exposure times of high resolution lUE 
spectra were rather long, i.e. ~ 2 to 3 hours, which is in fact 
comparable to the wind flushing time in these stars. We pro- 
vide here new perspectives on variability in the fast winds of 
CSPN by exploiting FUSE (AA905 to 1187 A) time-series 
observations. FUSE can uniquely deliver the requisite high 
signal-to-noise data in short integration times, while also 
potentially providing access to unsaturated resonance lines, 
including Pv, Siv, Svi, Ovi, and Cm. 



1.1 The target, NGC 6543 

We present in this paper a study of structure and variability 
in the fast wind of the central star of the planetary nebula 
NGC 6543 (HD 164963; BD +66° 1066). The nebula is well 
studied, with HST images in particular revealing complex 
structures that include phenomena such as bubbles and pre- 
cessing jets (e.g. Harrington & Borkowski 1994; Balick 2004; 
Wesson & Liu, 2004) . X-ray emission has been detected with 
Chandra (Chu et al. 2001; Guerrero et al. 2001), which is 
consistent with a point source at the central star and dif- 
fuse emission from within the elliptical central shell. Various 
studies have reported evidence for time-variable changes as- 
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Figure 1. Mean FUSE spectrum of the central star in NGC 6543 
from the F034 run in January 2007. 



sociated with the central star, including Mendez, Herrero 
& Manchado (1990; optical spectroscopy). Bell, PoUacco & 
Hilditch (1994; photometry), and Patriarchi & Perinotto 
(1997; UV spectroscopy). Fundamental parameters of the 
central star based on non-LTE model atmosphere analysis 
have been derived by de Koter et al. (1996; ISA- Wind code) 
and Georgiev et al. (2006; CMFGEN). Key adopted param- 
eters for the central star in NGC 6543 are listed in Table 1. 
We note that for a '/3=l'-type wind velocity law (e.g. Sect. 
5), the wind flushing time to travel from 0.1 to 0.9 of the 
terminal velocity (voo) is only ~ 45 mins. 



2 FUSE TIME-SERIES DATA 

Our study of the fast wind of NGC 6543 is largely based 
on FUSE data from Program F034 (P.I. - D.L. Massa), 
secured between 2007 January 13 to 16. A summary log of 
observations is provided in Table 2. The observations were 
obtained through the MDRS (4 x 20 ) aperture, spanning 
a total wavelength range of 905 to 1187 A at a spectral res- 
olution of ~ 15 km s"'^. (The instrument and performance 
is discussed by Moos et al. 2000 and Sahnow et al. 2000.) 
Typical individual integration times in HIST mode are ~ 



Clumped fast wind of NGC6543 3 



Table 2. Summary log of FUSE observations. 



Obs. ID 


No. of exp. 


UT Date 


MJD (start) 


AT (hours) 


F0340105 


26 


2007-01-13 


54113.1260 


9.7 


F0340106 


18 


2007-01-14 


54114.1097 


5.9 


F0340107 


19 


2007-01-15 


54115.0933 


7.3 


F0340108 


9 


2007-01-16 


54116.0680 


1.4 



10 minutes. Using the MDRS aperture unfortunately meant 
that occasionally the target was not fully captured and there 
were cases of light loss. The most usable components are 59 
spectra (processed through CalFUSE version 3.1.8) in the 
LiF2 channel (segment A; AA1086 to 1182 A) which covers 
the strategic line of Pv AA1117.98, 1128.01. A sparser se- 
quence of 29 spectra for the LiFl channel (segment A; AA987 
to 1082 A) provided some information on the behaviour of 
the (stronger) Ovi AA1031.92, 1037.62 wind-formed line. 

The mean FUSE spectrum is shown in Figure 1. Aside 
from Pv and Ovi, a strong P Cygni line is also seen in 
Svi AA933.38, 944.52. There is some evidence for blue-ward 
absorption in the excited Niv A955.34 lino profile and Cm 
A1175.67. The remainder of the spectrum is dominated by 
narrow interstellar and circumsystcm absorption lines due 
to atomic species and molecular hydrogen. 

Overall the mean FUSE spectrum of the central star 
in NGC 6543 is very similar to that seen in archival data 
in 2001 October (Prog. Q108) and March 2006 (prog. 
F034). The spectrum in Fig. 1 also closely matches the OR- 
PHEUS spectrum obtained in 1993 September by Zweigle 
et al. (1997). Note for completeness that the lUE high- 
resolution spectrum of NGC 6543 is dominated by satu- 
rated P Cygni profiles of CivAA1548.20, 1550.77 and Nv 
AA1238.82, 1242.80, plus well developed blueward absorp- 
tion in Ov A1371.29 and the excited line of Niv A1718.55 
(e.g. Patriarchi and Perinotto 1997). 



3 VARIABILITY 

The Pv resonance doublet in NGC 6543 is a well devel- 
oped but unsaturated P Cygni profile and thus offers the 
primary diagnostic of variable fast wind conditions in our 
dataset. A pair of Pv line profiles separated by ~ 2 days is 
shown in Fig. 2, to illustrate what is essentially the maxi- 
mum extent of profile fluctuations evident. A correction for 
a heliocentric systemic velocity of —66 km s~^ (Schneider et 
al. 1983) was applied by subtracting this value. Mostly the 
flux changes arc at ~ 10 to 20% of the continuum level and 
occur over localised blueward velocity regions, as opposed 
to the flux increasing or decreasing simultaneously over the 
entire absorption trough. As is the case with massive OB 
stars, the Pv P Cygni emission profiles in NGC 6543 are 
quite constant. The overall variance in the full 2007 time- 
series is certainly significant between ~ —500 to —1300 km 
s~^. The equivalent width of the Pv line is not conserved, 
and fluctuates between ~ 0.7 to 2.2 A (s.d. 0.3 A) for the 
blue component of the doublet (e.g. measured between 1111 
to 1120 A). 

Any changes in the other spectral lines are rather sub- 
tle if present at all. The lower panel in Fig. 2 shows the 
corresponding line profiles of Ovi A1037.62. Even the low 
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Figure 2. Fluctuations in the Pv resonance line profiles com- 
pared to the more stable low-velocity regions of Ovi A1037.6 



velocity (redward of ~ —600 km s^^) region of this line is 
likely optically thick and variance analysis does not support 
any changes in Ovi that are significant above a 95% confi- 
dence level. 

The data for Pv are by contrast very revealing. Rep- 
resentations of the variability evident are displayed in Fig- 
ure 3 as dynamic spectra (i.e. two-dimensional images of 
individual spectra ordered in time). The images reveal clear 
evidence for systematic and organised line profiles changes. 
Specifically, localised (in velocity) absorption enhancements 
are seen (darker shades in Fig. 3) migrating blueward from 
~ —400 to ~ —1000 km s~^. We identify at least 4 separ 
rate episodes of recurring features. The patterns in Fig. 3 
are essentially identical to the 'discrete absorption compo- 
nents' (DACs) commonly seen in UV resonance lines formed 
in the radiation-pressure-driven winds of massive OB stars 
(e.g. Kaper et al. 1996; Prinja, Massa & Fullerton, 2002). 
The lower left-hand panels in Fig. 3 focus on three of the 
more prominent DAC-like features in NGC 6543, with cen- 
tral velocities that increase between ~ —830 to —1010 km 
s~^ over 1.1 hours, ~ —510 to —1290 km s~^ over 6.0 hours, 
and ~ -700 to -1070 km s"^ over 1.3 hours. 

The DAC velocities above provide (linear) accelerations 
of between 3 x 10"^ km s"^ to 8 x 10"^ km s"^. These val- 
ues for NGC 6543 are up to a factor of 10 faster than typical 
acceleration rates measured for features in OB star wind 
(e.g. Prinja 1998). However the characteristic radial flow 
time of the wind (which scales as ~ R*/foo) in NGC 6543 is 
only ~ 5 minutes (Table 1), compared to ~ hours for O stars. 
The time-scales associated with the DACs in NGC 6543 are 
therefore significantly greater than the wind flushing time 
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Figure 3. Dynamic spectrum representations of organised, hourly variability in Pv AA1118,1128. In the upper plot the ordinate is the 
sequential spectrum number, and the panel on its right is a temporal plot showing the relative time of each exposure. The grayscales 
in the lower plot show the organised, hourly variability. Starting from the bottom, data in the pair of (lower) left-hand and right-hand 
panels span a total of ~ 1.22 days and 1.02 days, respectively. To enhance the contrast in those images the individual spectra have been 
normalised by the mean profile for the time-series. 



Clumped fast wind of NGC6543 5 




Figure 4. Power spectrum (in arbitrary units) for the Pv FUSE 
time-series. The arrow marks the primary peak corresponding to 
a period of 0.17 days in Pv (see Sect. 4). 

over the line formation region, which suggests that it is very 
unhkely the wind structures in the PN central star are due to 
processes entirely intrinsic to the fast wind. The DACs can- 
not, for example, be due to mass-conserved shells or blobs 
'riding' with the outflow from the star. As is widely believed 
for OB stars, the perturbations more likely arise as a re- 
sult of inhomogeneities close to the stellar surface, such as 
pulsations or magnetic fields. The empirical evidence pro- 
vided here suggests that the physical mechanism for initiat- 
ing wind structure in the fast wind of NGC 6543 may the 
same as that operating in massive, luminous stars (see Sect. 
6). 



4 MODULATED OR CYCLIC BEHAVIOUR 

In a few case studies based on very extended lUE spec- 
troscopy it has been demonstrated that the DACs in OB 
stars are modulated or quasi-cyclic, on variability time- 
scales that may relate to stellar rotation periods (e.g. Massa 
ct al. 1995, FuUcrton ct al. 1997, de Jong et al. 2001, Prinja 
et al, 2002). From the viewpoint of nebular shapes, cen- 
tral star binarity, and wind structure in line-driven winds, 
it is obviously interesting to consider whether the optical 
depth changes seen in the Pv lines of NGC 6543 are system- 
atically variable. The rotation period of NGC 6543 is not 
known, and determinations of the projected rotation veloc- 
ities are sparse for PN central stars in general. Note also 
that the central star in our study is a member of the ZZ Lep 
class of variable stars (e.g. Handler 2003). Bell, PoUacco & 
Hilditch (1994) document differences of ~ 0.01 mag in 5300 
A CCD photometry though with little evidence for periodic 
variations that are consistent over the four nights of their 
observations. Individual nights provided evidence for mod- 
ulations over ~2.6 hours and 3.2 hours. Analysing the same 
dataset, Handler (2003) cites a photometric modulation of 
3.5 ± 1.5 hours. The origin of the photometric changes in 
NGC 6543 may relate to photospheric (pulsational) changes 
or indeed the inner regions of the fast wind. The strong FUV 
photospheric lines in our FUSE data are relatively fixed in 
velocity (to within ± 10 km s~^) and do not support any 
significant role for binary motion. 

A periodogram analysis was carried out to search for ev- 
idence of repetitive or cyclic properties in the Pv line profile 
changes. The basic Fourier method involved uses the iter- 
ative CLEAN algorithm (Roberts, Lehar & Dreher 1987) 
to deconvolve the features of the window function from the 
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Figure 5. Individual Pv spectra are shown phased over two cycles 
on a period of 0.17 days; a clearly coherent behaviour is evident. 
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Figure 6. Variance spectrum for low-resolution lUE time-series 
of the Civ AA1548,1550 line in NGC 6543. The horizontal dashed 

line represents the 95% confidence limit. The lower panel shows 
the low-resolution line profile (solid line) together with the a rep- 
resentative high-resolution lUE spectrum (dashed line). 
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Figure 7. SEI model fits to the mean Svi AA933, 945 and Pv AA1118,1128. The adopted optical depth bins as a function of normalised 
velocity (= v/voo) are shown in the right-hand panels. 



discrete Fourier transform. (A gain of 0.5 with 100 itera- 
tions was used.) The power spectrum sampled over ~ —600 
to —1100 km (i.e. the range of well-developed DACs) in 
Pv is shown in Fig. 4. The data sampling is generally not in- 
tensive enough for a very robust period search and the power 
in Pv is spread across several frequencies. Nevertheless we 
identify as potentially interesting the main power peak in 
Pv at ~ 5.85 days~^, corresponding to a period of ~ 0.17 
day. The half-width at half-maximum of the zeroth order 
peak in the window function provides an approximate (and 
conservative) estimate of the error in the frequency of ~ 0.1 
day~^. Independent determinations of the period sampled 
over 300 km s~^ bins in the region between —600 to —1100 
km s~^ suggest an error in the frequency of ~ 0.5 day~^. 
The adjacent peaks in the power spectrum (Fig. 4) are com- 
paratively less stable in similar velocity bin sampling. 

Grayscales of individual spectra normalised to the mean 
arc shown in Fig. 5 phased on the 0.17 day period. The 
diagnostic Pv data do clearly show some coherency on the 
0.17 day time-scale. The modulation is represented by the 
occurrence of two sequential episodes of migrating structures 



in the wind. (Recall that the total length of the FUSE time- 
scries combined here is ~ 3 days, and therefore spans several 
'cycles'.) 

Tentative evidence of variability on hourly time-scales 
can also been gleaned by examining low-resolution (R ~ 300) 
lUE spectra of the central star. We have retrieved from 
the ESA lUE archive a sequence of 38 spectra (SWP54865 
to SWP54960), which extend over 4.2 days. Despite being 
saturated, the blue-wing of Civ is apparently highly vari- 
able (e.g. Patriarchi & Perinotto, 1997). Furthermore, stud- 
ies of fluctuations in the extended blue-wings of saturated 
P Cygni profiles have (perhaps surprisingly) revealed evi- 
dence for modulated behaviour (see e.g. Howarth, Prinja & 
Massa, 1995; Prinja et al., 1998). We used high-resolution 
lU E spectra of NGC 6543 that were degraded (by convolu- 
tion with a 5A full- width at half-maximum Gaussian), to im- 
prove the internal consistency of the low-resolution data and 
to normalise the continuum fluxes near Civ. The variance 
characteristics of the low-resolution lU E Civ time-series is 
shown in Fig. 6. There is a double-peaked nature to the 
TVS (e.g. Fullerton, Gies & Bolton 1996) spectrum, and 
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we elect to sample in particular the profile changes between 
— 1000 and —2000 km s~^, since these are likely to represent 
variability in the extreme blue-wings. However, no signal is 
detected in CiV at the frequency identified for Pv modula- 
tions. The lUE data quality prevent us from making any 
firm conclusions £is to whether any (quasi-)periodic signal 
in NGC 6543 is long-lived. We also attempted to phase-fold 
the original photometric data of Bell et al. (1994, obtained 
in 1993; Don Pollacco private communication) but did not 
find any convincing evidence for a modulation on ~ 0.17 
days. 



5 LINE-SYNTHESIS AND OPTICAL DEPTHS 

In order to further quantify the physical properties of the 
DACs in the wind of NGC 6543, the FUSE resonance Unes 
were modelled using the SET method (see Lamers, Cerrutti- 

Sola & Porinotto 1987), with the modified procedures of 
Mcissa et al. (2003). The radial optical depth, rrad(ti'), is 
treated as 21 independent velocity bins (each ~ 0.05 Voo 
wide) to provide greater flexibility in matching absorption 
optical depths in a structured wind that is affected by vari- 
able features such as DACs. We assume a standard param- 
eterised velocity law of the form. 
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(1) 



where w = v/voo and x = r/Ri,. The effect of a photo- 
spheric spectrum in the SEI profile fitting is approximated 
by a Gaussian of full- width at half-maximum — 250 km s~^; 
its effect is only to improve the match at low velocities (^ 

0. 2 Voo) and does not impact seriously on the overall wind 
parameters derived in this section. Model fits to the mean 
Pv profile in our time-series is shown in Fig. 7, where a very 
good overall match is achieved for Voo = 1400 km s^^, ,5 = 

1, and Wturb (the small-scale velocity dispersion parameter) 
= 0.05. Note that the redward component of Pv may be af- 
fected by the (neglected) contribution due to Siiv AA1122.5, 
1128.3. These values are also consistent for reproducing the 
blue- wing of a much stronger (saturated) profile such as Svi 
(Fig. 7). 

The central star is H-rich, with normal solar abundance 
ratios though Carbon may be overabundant (e.g. Georgiev 
et al. 2006). The SEI model of Pv in Fig. 7 therefore gives a 
(mean) value for the product of mass-loss rate and ionization 
fraction, M q{P'^^), integrated over 0.2 to 0.9 Voo,o{^ 5.0 x 
10~^ Mq yr~^ (for parameters in Table 1). Over the velocity 
range where the peak variance occurs in Pv (i.e. ^ 0.5 to 
0.9 Voo), SEI fits to the cases of minimum and maximum 
observed absorption yield M q{P'^^) values of ~ ~ 1.9 x 
lO"'' Mq yr-i and 3.5 x lO"'' Mq yi'^. 

The fact that the SEI model can match well the over- 
all morphology of the observed absorption in Pv permits 
us to model individual line profiles in the FUSE time- 
series. Specifically we focus on the data covering the two 
most prominent DAC-like features in lower left-hand panels 
of Fig. 3 (in both cases the features start from the nor- 
malised time of 0.0 hours; the corresponding FUSE expo- 
sures are F03401050012 to F03401050082 and F03401060042 
to F03401060122). The fits to these spectra, which are vary- 
ing due to the occurrence of DACs, can then be compared 
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Figure 8. SEI model fits (dashed lines) to sequences of Pv line 
profiles where two of the prominent DACs are migrating between 
~ —600 km to —1100 km s~^; i.e. the lower left-hand panels 
in Fig. 3. The isolated profile fitted in the upper left panel is a 
representative 'no-DAC (maximum flux) case. 
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Figure 9. Ratio of optical depth in the DAC to minimum optical 
depth at that velocity from SEI fits, for two independent episodes 
(i.e. open and closed symbols) identified in Figs. 3 and 8. 



with a 'no-DAC case i.e. a spectrum representative of max- 
imum flux over most velocities (i.e. F03401080072). 

The SEI models fits to the Pv (time-series) lino profiles 
are shown in Figure 8. The fiexible optical depth bins are 
designed to allow matches to a wind profile that is obvi- 
ously not 'smooth', but is instead distorted due to density, 
ionization or velocity features (see e.g. Massa et al. 2003). 
The good overall reproduction in Fig. 8 permits us to ex- 
tract optical depths at the specific velocity of the migrating 
DACs, and normalise them to the no-DAC case. Ratios of 
T(DAC)/T(min) (obtained from fits to the individual pro- 
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Figure 10. Key diagnostics lines in the FUSE and high-resolution lUE spectra (dots) of the central star are shown with CMFGEN 
model fits for /oo = 0.06 (solid line), 0.08 (dashed line) and 0.1 (dotted line). Notice the particular sensitivity of Pv to the clumping 
factor. 



files) arc plotted as function of normalised velocity (v/voo) 
in Fig. 9. The plot therefore shows the optical depth evolu- 
tion of two independent DAC episodes as a function of their 
central velocities. In both cases the optical depth ratio due 
to DACs rises to a maximum of ~ 5 at ~ 0.62 Vao , and drops 
to ~ 2 toward 0.75 Voo. The indication that the maximum 
contrast of the DACs is reached at intermediate velocities 
is somewhat unexpected. The increase in Sobolev optical 
depth at these positions may bo duo to a plateau in the wind 
velocity law, i.e. a radially extended region of near-constant 
velocity. Wind material flows through the plateau, and the 
observed accelerations of DACs are then associated with the 
evolution (in velocity and radial extent) of the plateau as a 
function of time. The origin of such a plateau could, for 
example, relate to the formation of co-rotating interactions 
regions in the outflow (see e.g. Cranmer & Owocki, 1996, 
and discussion in Sect. 6). 

To gain some insights into the ionization balance and 
wind clumping in NGC 6543, we also employed the unified 
model atmosphere code CMFGEN (Hiller & Miller, 1998), 
which calculates non-LTE models of spherically symmetric 
winds with multi-level atoms. The fundamental parameters 
listed in Table 1 were adopted, together with log g = 4.2 and 
a He/H abundance ratio (by number) of 0.1 (e.g. Georgiev 
et al. 2006). A parametric treatment of wind clumping is 



implemented in CMFGEN, in terms of the volume filling 

factor, /, defined as 

/ = /oo + (l-/oo)ea;p(-«KO, (2) 

where Vd is the velocity at which clumping starts, and which 
we adopt as 35 km s~^ (i.e. above the sonic point). Our ob- 
jective here is primarily to investigate the spectral signatures 
of clumping in CMFGEN models applied to the fast wind of 
NGC6543; we do not attempt a detailed exercise in fitting 
several lines. 

Models were run for cases of the clump volume filling 

factor, /oo ~ 0.06, 0.08 and 0.1, and the synthesised spectra 
are compared in Fig. 10 to key wind lines in the FU SE and 
high-resolution lUE data of the central star. For a mass- 
loss rate of 6 x 10~* M0 yr~^ the models provide good 
overall fits to the observed resonance and excited state line 
profiles. It is interesting to note that the Pv lines arc par- 
ticularly sensitive to the clumping factor, more so than Ov 
and Niv. Our data favour the case of /oo = 0.08, though 
of course we neglect here the consequences of these param- 
eters for reproducing diagnostic optical lines in the central 
star of NGC 6543, such as Hell A4686 and Civ A5801. Figure 
11 displays the predicted ionization fractions of P^+, P''"'", 
and P®"*" as a function of velocity, and for the three different 
clumping factors. The dominance of P^"*" at low velocities re- 
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Table 3. Comparison of NGC 6543 and mid-O star DAC properties. 



DAC measure 


NGC 6543 


O stars (e.g 


07 III) 


^initial 


~ 0.3 Voo 


^ 0.3 Doc 




Recurrence 


yes (over ~ hours) 


yes (over ~ 


days) 


{dv/dt)x{R*/Doo) 


~ 12-24 km s-l 


~ 10 km s- 


1 


maximum r(DAC)/r(min) 


~ 5 


~ 5 




Quasi-periodic 


yes; P ~ 0.14 days 


yes; e.g. P " 


-J days 


FWHM/'Uoo (initial) 


0.15 


~ 0.3 Doo 




FWHM/doo (final) 


^ 0.05 Voo 


< 0.1 Voo 





suits in the mis-match with the observed profile seen in this 
region (Fig. 10). The ion fractions of and P*"'" increase 
over all velocities for lower values of foo, with an increased 
recombination from P^"*" to P^"*" particularly evident at the 
highest velocities. The CMFGEN models run here generally 
predict C*^ and P'''^ as dominant (i.e. qi > 0.97 over ~ 0.2 
to 0.9 Voo), and and C"*" close to dominant (i.e. qi ~ 
0.6 and 0.2, respectively). 

Our results indicate that Pv is a sensitive spectral di- 
agnostic of small-scale clumping in the fast wind, together 
with the large-scale structures represented by the DACs dis- 
cussed in Sect. 3. We highlight the need to adopt fairly small 
volume filling factors in order to match the Pv resonance 
line doublet. P*+ is a trace ion, just below dominant (Fig. 
11), and is therefore extremely sensitive to changes as pa- 
rameterised here, since its abundance depends strongly on 
density. This is in contrast to the results for O stars pre- 
sented by Bouret et al. (2005), where P*"*" is dominant and 
not as sensitive to small fluctuations in density. 



6 DISCUSSION 

The empirical evidence presented here suggests that (i) the 
fast wind of NGC 6543 is highly variable on time-scales 
of hours, (ii) the wind changes are systematic in terms of 
recurrent absorptive optical depth structures that migrate 
blueward, and (iii) the temporal behaviour of the Pv res- 
onance line P Cygni profiles is very similar to that due to 
the presence of 'discrete absorption components' (DACs) in 
luminous OB stars. A comparison between the DAC prop- 
erties of NGC 6543 and a 'representative' mid-0 giant is 
presented in Table 3, where the latter is summarised from 
various studies (see e.g. those cited in Sect. 3). We note that 
the fundamental characteristics such as initial detection at 
low Doppler velocities, recurrence, modulation, and optical 
depth contrast are broadly similar for both cases; the same 
is also true for the (linear) acceleration scaled by the flow 
time-scale (R*/uoo). There is some indication that the line- 
of-sight velocity dispersion of the DACs in NGC 6543 is 
somewhat smaller than for a typical 0-type giant (such as 
^ Per or 68 Cygni). 

Our interpretation of these results therefore is that the 
physical mechanism responsible for the DAC structures in 
the fast wind of NGC 6543 is likely the same as that acting 
in massive OB stars. There remain several uncertainties in 
our understanding of this phenomena in massive stars, but 
in the particular case where quasi-periodic wind variabil- 
ity is detected, several authors have promoted a scenario in 
terms of co-rotating interaction regions (CIRs; e.g. Cran- 




500 1000 

Blueward velocity (km s ^) 

Figure 11. Log of the PiV, Pv and Pvi ion fractions as a function 
of blueward velocity, for cases of foo = 0.06 (solid line), 0.08 
(dashed line) and 0.1 (dotted line). 



mer & Owocki 1996; FuUcrton et al. 1997; dc Jong et al. 
2001). The notion is that the DACs are spectroscopic signa- 
tures of spatial structures in the wind that are causally con- 
nected to stellar surface irregularities. These photospheric 
inhomogeneities cause the wind from different longitudinal 
sectors on the stellar surface to emerge with different densi- 
ties and/or velocities. The consequence is to form different 
adjacent streams that meet to create spiral-shaped CIRs. 
Key observational properties of the DACs can be matched 
with this model, with variable optical depth enhancements 
arising from the combination of a plateau in the radial veloc- 
ity as well as a density perturbation. In this model, the vari- 
ations in absorption which occurs over significantly larger 
time-scales than the wind fiow speed are due to gas fiow- 
ing through a semi-permanent, co-rotating structure. The 
rotation rate of the central star in NGC 6543 is unfortu- 
nately not known, and we are unable to comment further 
in this context on the significance of the ~ 0.17-day pe- 
riod identified in Sect. 4. In the case of OB stars the re- 
currence of DACs is mostly only quasi-cyclic, and has only 
rarely been demonstrated to be an integer relation to the 
stellar rotation rate. Nevertheless, our results provide the 
constraint that the mechanism for forming coherent pertur- 
bations in the outflows is apparently operating equally in the 
radiation-pressure-driven winds of widely differing momenta 
{MvooR*'^) and flow times, as represented by OB stars and 
CSPN. 
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As regards the 'trigger' provided by photospheric ir- 
regularities for the growth of wind structure, the action of 
surface velocity fields due to pulsation and/or ordered mag- 
netic fields may bo just as applicable to CSPN as it is to 
OB stars. Multi-period pulsation modes arc thought to bo 
the cause of short-period (~ hourly) variability seen in the 
photometric measurements of several central stars, including 
both H-rich and H-deficient cases (e.g. Handler et al. 1997; 
Ciardullo & Bond 1996) . Kilogauss strength magnetic fields 
have also been detected (using spoctropolarimctry) that are 
associated with the PN central stars (e.g. Jordan, Werner 
& O'Toole 2005), where they may be located mostly in the 
envelope, as opposed to the degenerate core, and can be 
affected by mass-loss. One caveat here is that the link be- 
tween pulsation or magnetic fields and repetitive wind struc- 
tures in OB stars has thus far proved observationally very 
challenging and generally unconstrained; establishing such 
a connection in CSPN is going to be even more difficult! 

The ab inito CMFGEN calculations presented in Sect. 
5 favour a small clump volume filling factor, foo ~ 0.08. 
Whether or not the parameterisation of clumping incorpo- 
rated in CMFGEN is physically meaningful, we note that 
the high sensitivity of the Pv lines to the clumping factor is 
again similar to that apparent for the same line in O-type 
stars. Following the clumping studies of e.g. Markova et al. 
(2004), Bouret et al. (2005), FuUerton et al. (2006), our re- 
sults for NGC 6543 would imply a reduction by a factor ~ 4 
(i.e. 1 / V foe ) in the mass- loss rate duo to the clumped fast 
wind compared to smooth-wind models. The analysis of op- 
tical lines with the non-LTE code FASTWIND by Kudrtizki, 
Urbaneja & Puis 2006) also points to substantially clumped 
CSPN fast winds. Downward revisions in the central star 
mass-loss rate would affect the dynamics of the hot central 
cavity excavated by the feist wind, thus modifying its interac- 
tion with the surrounding nebula, including contributions to 
the X-ray emission. Naturally our far-UV analyses needs to 
be extended to cover a large sample of central stars to permit 
firm assessments of any revisions to CSPN mass-loss rates, 
and most importantly, tested for consistency between UV 
and optical diagnostics of clumping. Finally, looking ahead, 
the vastly increased sensitivity soon to be offered at radio 
wave-bands by the upgraded e-MERLIN and EVLA facilities 
raises new possibilities for measuring the free-free emission 
from the fast wind itself and thus provide a powerful new 
constraint on the central star mass-loss. 
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